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Abstract 


The problem of the cosmic coincidence is a longstanding puzzle. This 
conundrum may be solved by introducing a coupling between the two dark 
sectors. In this Letter, we study a coupled quintessence scenario in which 
the scalar field evolves in a power law potential and the mass of dark matter 
particles depends on a power law function of cf- It is shown that this scenario 
has a stable attractor solution and can thus provide a natural solution to the 
cosmic coincidence problem. 
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Before the accelerated expansion of the nniverse was revealed by high red-shift snpernovae 
la (SNe la) observations [1], it conld hardly be presnmed that the main ingredients of the 
nniverse are dark sectors. The concept of dark energy was proposed for nnderstanding 
this cnrrently accelerating expansion of the nniverse, and then its existence was confirmed 
by several high precision observational experiments [2,3], esp. the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite experiment [2], The WMAP shows that dark energy 
(DE) occnpies abont 73% of the energy of onr nniverse, and dark matter (DM) abont 23%. 
The usnal baryon matter which can be described by onr known particle theory occupies only 
about 4% of the total energy of the universe. Although we know that the ultimate fate of the 
universe is determined by the nature of DE, the information about its nature we can acquire 
is still very limited. So far the conhrmed information about DE can be summarized as the 
following three items: it is a kind of exotic component with negative pressure such that it 
can accelerate the universe’s expansion; it is spatially homogeneous and non-clustering; and 
it dominates the universe today although it contributes little to the universe at the early 
times. 

The investigation of the nature of DE is an important mission in the modern cosmology. 
Much work has been done on this issue, and there is still a long way to go. One candidate for 
DE is vacuum energy density or cosmological constant A [4] for which the equation of state 
w = —1. The cosmological model that consists of a mixture of vacuum energy and cold dark 
matter (CDM) is called LCDM (or ACDM). Another possibility is the so-called QCDM 
cosmology which based upon a mixture of CDM and quintessence held [5]. The energy 
density and the negative pressure are provided by the quintessence scalar held 0 slowly 
evolving down its potential E(0). The equation of state of the quintessence —1 < tc < —1/3 
is guaranteed by the slow evolution. However, as is well known, there are two difficulties arise 
from all of these scenarios, namely, the “hne-tuning” problem, and the “cosmic coincidence” 
problem [6]. The cosmic coincidence problem states: Since the energy densities of DE 
and DM scale so diherently during the expansion of the universe, why are they nearly 
equal today? To get this coincidence, it appears that their ratio must be set to a specihc, 
inhnitesimal value in the very early universe. 

A possible solution to this cosmic coincidence problem may be provided by introducing a 
coupling between quintessence DE and CDM. This coupling is often described by the varible- 
mass particle (VAMP) scenario [7]. The VAMP scenario assumes that the CDM particles 
interact with the scalar DE held resulting in a time-dependent mass, i.e. the mass of the 
CDM particles evolves according to some function of the scalar held 0. It has been shown 
that if we choose an exponential potential to the quintessence scalar held 0 and at the same 
time assume the CDM particle mass also depends exponentially on 0, the late-time ratio 
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between DM energy density and DE energy density will remain constant [8]. This 
behavior relies on the existence of an attractor solution, which makes the effective equation 
of state of DE mimic the effective equation of state of DM at the late times so that the late 
time cosmology insensitive to the initial conditions for DE and DM. In addition, the hnal 
effective equation of state of both components is negative and may lead, if this value is less 
than —1/3, to an accelerated expansion of the universe. Therefore, the scenario containing 
quintessence with exponential potential and VAMPs with exponential function of 0 solves 
the cosmic coincidence problem in this sense. However, it should be pointed out that this 
case is not the unique case. We consider here another case — the power law case. 

In this Letter, we explore the case in which the scalar held evolves in a power law 
potential and the mass of VAMPs depends on a power law function of 0. We hnd that the 
power law case also has a stable attractor solution which is similar to the exponential case. 
We numerically solve the equation of motion of the scalar held 0 and then illustrate the 
behavior of DE and DM. The numerical results show that the power law case of coupled DE 
with VAMPs also solves the cosmic coincidence problem. 

Consider, now, the CDM particle y with mass M depending on a power law function of 
the DE held 0, 


M^(0) = , (1) 

where 0 is expressed in units of the reduced Planck mass Mp {Mp = I/vottG = 2.436 x 
lO^^GeV), and a is a positive constant. The scalar held has a power law potential 

^( 0 ) = Vy , ( 2 ) 

where /? is a positive constant. Since the CDM particle is stable, its number density must 
obey the equation 


= 0 , (3) 

where the dot denotes a derivative with respect to time, H = a/a represents the Hubble 
parameter, and a{t) is the scale factor of the universe. The energy density of DM p^ is also 
0-dependent, which is given by 


p^(0) = , 


and it follows that 


Px A ^^Px ~ ^ A Px- 


( 4 ) 

( 5 ) 
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Since the total energy of DE and DM is conserved, the equation of motion for DE can be 
obtained 


+w^) 



( 6 ) 


where the usual parameter of equation of state for the homogeneous scalar held is given by 

^ ^ ^ -v{4>) 

p^~ + V{(j>) ■ 

The equations of motion (5) and (6) can also be written in the form of effective equations 
of state for DM and DE 


where 


Px + 3^Px(l + = 0 > 

+ 3Hp^{l + = 0 , 


,(e) _ 

^ 3H(t) 


a (j)' 
3^ ’ 


( 8 ) 

(9) 


( 10 ) 



a ^Px 
3H (pp^ 


^ 0 Px 

''-'*-3^7, 




( 11 ) 


are the effective equation of state parameters for DM and DE, respectively. Primes denote 
derivatives with respect to m = ln(a/ao) = — ln(l + z), where 2 : is the red-shift, and Qq 
represents the current scale factor. From (6) or (9) one can get the equation of motion for 
the scalar held 0, 


0 -1- 3Hp 


a 


-Px 


.tv 


( 12 ) 


The Friedmann equation for a spatially hat universe with DE, DM, baryons, and radiation 
reads 

3H^ = —0^ -|- V (0) + Px{4>) + Pb + Prad , (13) 

the Planck normalization Mp = 1 has been used here. Using the Friedmann equation (13), 
the equation of motion of held p can be written in terms of the derivatives of p with respect 
to u, 

+ ^(Px + Pb + gPrad + 2U)0' = —p^ — , (14) 
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where 


j^2 _ 1 (Px + Pb + Prad + V) 

“3 1-0'V6 


(15) 



FIG. 1. The phase plane for a = 11 and /3 = 4. The units of K and M* are in pco and 
Pcoln^o-, respectively. The three lines correspond to cases (14, M*) taken to be (0.1,230), (0.2,23) 
and (0.3,2.3), respectively. 


It is known that the attractor solution can be given analytically in the exponential case. 
In this case, however, it is difficult to obtain the analytical solution of the attractor. The 
reason arises from that the late-time behavior of the held 0, unlike in the exponential case, 
is not linear in this power law case. In what follows we will show this difficulty. At the 
late times, the energy densities of baryons and radiation are red-shifted away completely, 
there remains only DE component and DM component if the attractor exists. This requires 
in this limit, it follows that 





^ i/2 ^ 0 

3*^ ~ 3^ ■ 


(16) 


The attractor makes the ratio between energy densities of DE and DM be a constant, i.e. 
= 0, which gives 


202- a00" + a0'2 = 0 . 


(17) 


We see explicitly that the situation in the case of power law is more complicated than in the 
case of exponential. In the exponential case, 0' = const., then using the Eriedmann equation 
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this constant can be determined, furthermore it can be shown that and are both 0 - 
independent constant.^ However, in this case, since 0' 7 ^ const., we can not identify the 
attractor easily. We are now required to solve the equation of motion of 0 (14) numerically. 
Whether the attractor solution exists or not can be conhrmed by the numerical results. 



FIG. 2. Top panel: A typical solution for the differential equation (14), namely the evolution 
of the 0 field. The corresponding parameter configuration is; a = 11, 0 = 4, 10 = O.lpcO) and 
M* = 230/9co/n-xO- Bottom Panel: the evolution of 0' for the same parameters used in top panel. 


The numerical results indeed show that there exists a stable attractor solution which 
depends only on the parameters a and 0 while is very insensitive to the initial conditions 
and the chosen values of 10 and M*. As an example, we take the case o; = 11 and 0 = 4. 
The (0,00 phase diagram is plotted in Fig.l. It is shown in the phase plane that lines 
corresponding to different chosen values of ( 10 , M*) will converge together to the attractor 
solution with the cosmological evolution. Note that we express 10 and M* in units of pco and 
Pco/n^Oi respectively, where pco = 4.2 x lO^'^^GeW^ is the present critical density of universe. 
It is of interest to notice that the lines in the phase plane will undergo a period of oscillation 
before they converge to the attractor. From Fig.l, we see that in the attractor region 0 and 


^As an example, we take M^(0) = M*e , and F(0) = then the attractor solution is 

given by 0 = 00 + = 1 - It can be seen that the attractor depends only 

on A and g. 
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0' both tend to 0 and at the same time 0V0 becomes a constant. Using the fact that in 
the attractor region dVeff/dcj) ~ 0, where V^ff = U + is the effective potential, one gets 
ap^ ~ PV in this limit. We re-express this condition as 

1 


Considering hly = 1 — and 


Combining with (16), we get 


0 , one obtains 
a 




Ct -]r (d 


(f) a + P ^ 

consequently the attractor solution in the held space can be expressed clearly as 


— 0ne 


y+l3 


such that 


U'W = <> = 


a 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


cy P 

The existence of the attractor solution indicates that the role the coupled quintessence model 
in a power law case plays in solving the cosmic coincidence problem is similar to that of the 
model in an exponential case. 
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FIG. 3. Top panel; The evolution of the relative abundance of different species, expressed 
as fractions of the critical density. The corresponding model parameters are: a = 11, /3 = 4, 
K = O.lpcO) and M* = 2‘i{)pco/n^o. Bottom panel; Effective equations of state for DE (solid line) 
and DM (dashed line) for the same parameters used in top panel. 

We now illustrate the cosmological consequence of this scenario. The solution of the 
differential equation (14), namely the evolution of 0 (m), is plotted in the top panel of 
Fig.2. Note that for concreteness the values of (14, Tf*) are taken as 14 = O.lpco and 
M* = 230pco/''^xo- explicitly in this case that 4 increases in the beginning and begins 

to decrease around u ~ —2.5, which implies that the 0 held climbs slowly up the potential 
initially and then rolls down with the evolution of the universe. Moreover, we notice that 
around m ~ 0 some wiggle appears. To make the picture of the kinematics of the held more 
clear, we plot the Velocity’ of the held 4 '(m) (note that 0' = <p/H) versus u = — ln(l + z) 
in the bottom panel of Fig.2. It is shown that the wiggle in the 0-evolution diagram corre¬ 
sponds to the rapid oscillation of the rolling-down velocity. After the period of oscillation 
the held system will enter a stable attractor regime. 

It is remarkable that the feature of the coupled quintessence model is quite diherent 
from the scalar DE model in the absence of the interaction with DM. In Fig.3, we show the 
evolution of density parameters of various components in universe and the ehective equations 
of state for DE and DM. The model parameters are taken to be the same as above. It is 
exhibited that after a transient period of baryon domination, the universe enters into a 
DE-dominated epoch, and the interaction between DE and DM forces the ratio between 
energy densities of DE and DM to remain constant via DM particles varying mass. In the 
future, the universe will be governed by DE and DM which are coupled together, and pb 
and prad are diluted away with the usual laws, a~^ and respectively. We see that in 
the attractor regime the effective equations of state for DE and DM converge together and 
both are negative. Notice that the attractor solution is going to be reached today in this 
example. 

In summary, we investigate in this Letter the attractor solution of a coupled quintessence 
scenario in which the scalar held evolves in a power law potential and the mass of VAMPs 
depends on a power law function of 0. It has been shown by numerical calculation that 
this interacting scalar held scenario has a stable attractor and can thus provide a natural 
solution to the cosmic coincidence problem. 
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